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ABSTRACT 

Aims. We present a detailed mass reconstruction of the cluster RCS0224-0002 at z = 0.773 from the strong lensing 
features observed with HST/WFPC2. 

Methods. The mass profile is reconstructed using a parametric approach. We introduce a novel method to fit extended 
multiple images based on the Modified Ifausdorfi^ Distance between observed arcs and the arcs reproduced by the model. 
We perform the detailed error analysis of the model parameter using the MCMC method. 

Results. Our model reproduces all the observed strong lensing features of the RCS0224-0002 and predicts the redshift of 
one of the arcs systems to be z ~ 2.65 (the other system has an spectroscopic redshift of 2; = 4.87). The reconstructed 
inner mass profile is well fitted by a non-singular isothermal sphere, rather than with an NFW model. Dark matter 
substructure, derived from the light distribution of the most luminous cluster members, is crucial for reproducing the 
complexity of the quadrupole image system, which could not be achieved otherwise. The reconstructed mass distribution 
closely follows the light, however it is significantly shifted from the X-ray emission of the gas. The mass of RCS0224- 
0002 derived from the lensing model, « 2 x lO^** Mq is in a very good agreement with the one obtained from the X-ray 
temperature measured with deep Chandra observations. 
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1. Introduction 

The presence of dark matter is evident in galaxy clusters, 
where the mass to light ratios, M/L « 200 h (Mq/Lq), 
by far exceeds the mass of stars, an d the hot gas implied 
by X-ray data ([Hradeckv et al.l[2000h . Gravitational lensing 
is one of the most attractive methods to directly study the 
mass distribution in the Universe on different scales, regard- 
less of its component and dynamica l state (see e.g. ESA- 
ESO Working Groups, Report No. 3, iPeacock fc Schneided 
12006 ^. Strong gravitational lensing in galaxy clusters leads 
to the for mation of multiple images and giant luminou s 
arcs (e.g. [Schneider. P. and Ehlers J. and Falco E.I I1992D . 
Observations of these features allow us to investigate the 
distribution of the mass responsible for the deflection, and 
in particular provide accurate estimates of the total mass 
within giant arcs and the innermost density profile, which 
can be compared with N-body simulations. In addition, 
mass models can reveal the location of highly magnified 
background galaxies and allow spectroscopic studies of very 
distant faint, highly magnified sources, that would be un- 
der normal circumstances beyond the reach of lOm-class 
telescopes. 

Accurate estimates of the mass profiles of galaxy clus- 
ters are fundamental for the modern cosmology since they 
provide severe tests for the theories of structure formation. 
The most popular theory based on the assumption of non 
interac ting cold dark matt er predicts a universal profile 
(nfw,|n avarro e t al. 1996) rather then an isotherm al pro- 
file. For the lensing cluster Abell 1689, ,Broadhurst et al.l 
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(|2005a[ ) ruled out the isothermal profile with lOcr con- 
fidence. However, a nother study of the same cluster by 
iHalkola et all (|2006[ ) shows that both the elliptical NFW 
and the isothermal softened elliptical fit the data well. 
Hence, although deep observations of this cluster showed 
a formidable arc system, parametric strong lensing models 
surprisingly lead to different mass density profiles, high- 
lighting the difficulty of current inversion techniques in de- 
termining the uniqueness of the solution and the real un- 
certainties of the reconstructed mass maps 

Here we present a study of the cluster RCS0224-0002 
at z = 0.773 which was discover ed as a part of th e Red- 
Sequence Cluster Survey (RCS. Gladders et al.ll2002D . After 
the identification of the main strong lensing features of this 
cluster with VLT spectroscopy, fo llow-up observa. t ions w ere 
carried out with HST-WFPC2 bv lGladders et all (l2002fl. i n 
X-rays with the Chandra observatory ( Hicks et al. 200qh 
and in sub-mm using SCUBA on the JCMT fWebb et all 
I2005D . 

We construct a parametric model of the projected mass 
density distribution of RCS0224-0002 based on its strong 
lensing features, one of which with secure redshift. The 
method used in this paper to construct the best mass 
model i s based on the so-called Modified Hausdorff Distance 
(MHD. " Pubuisson fc Jainl '1994'). and has the advantage of 
allowing us to use the information provided by the sub- 
arcsecond morphology of arcs. We compare the mass distri- 
bution with the spatial distribution of the hot gas obtained 
form the X-ray data. 

When we were finalizing this paper, a lensing model 
of the same cluster has been independently presented by 



2 



J. Rzepecki et al.: Strong Lensing Analysis of the Cluster RCS0224-0002 at z = 0.77 




ISwinbank et al.l (|2007f) . However, these authors focus their 
work on the properties of a highly magnified z = 4.87 
galaxy observed in the field; moreover, their lensing model, 
which is based only on the constraints provided by a single 
arc system (the giant arc labeled A in Fig. [T|), is signifi- 
cantly different from ours. 

The paper is structured as follows: In Sect. [2] we present 
data available on the RCS0224-0002. In Sect. [3] we discuss 
the strong lensing features and the red sequence of the 
RCS0224-0002. Section g] is dedicated to the X-ray emis- 
sion of the RCS0224-0002. In Sect. [S] we present assump- 
tions behind our model and the method we use to obtain 
the projected mass distribution. In Sect. [S]we present and 
discuss our results. In Sect. [7] we perform error estimation. 
And finally in Sect. [8] we present our conclusions. 

In this paper we use a standard cosmological model with 
r2,„ = 0.3, 17a = 0.7, and Ho = 72 km s'^ Mpc'^ We 
give all the magnitudes in the AB system, if not otherwise 
specified. 



2. Observations 

The HST observations of the RCS0224-0002 were taken on 
the 2001/08/20 in two filters, F606W and F814W using 
the WFPC2 camera (PI: Gladders, Proposal ID: 9135). The 
target coordinates were RA: 02:24:30.82, DEC: -00:02:27.8 
and the exposure time for each filter was 1100 seconds. 
The WFPC2 data reduction was performed by Associations 
Science Products PipelineQ 

The X-ray data were taken on the 2002/11/15 with 
the ACIS-S instrument on the Chandra observatory (PI: 
Gladders, Proposal Num: 03800013). The target coordi- 
nates were RA: 02:24:34.10, DEC: -00:02:30.90 and the 
exposure time was 14560 seconds. On the 2004/12/09, 
RCS0224-0002 was observed with the ACIS-S again (PI: 
EUingson, Proposal ID: 05800899) with exposure time of 
90150 seconds. The two ACIS-S observations were com- 
bined with CIAO 3.3, using CALDB 3.2.1, leading to 100.8 
ksec of effective exposure time. Details on the reduction 
and spectra l analysis, whose res ults are given below, can 
be found in iBalestra et al.l (|2007f) . 



3. Arc identification and cluster members 

RCS0224-0002 has seven prominent luminous arcs and ar- 
clets marked as Al, A2, A3, Bl, B2, B3, and B4 in Fig. [T] 
Unfortunately, out of those seven arcs, only one arc sys- 
tem (A) has a confir med spectroscopic redshift of 4.87 
(iGladders et al.l I2002D . The same authors estimated the 
redshift of system B within the range 1.4 to 2.7 based on the 
lack of emission lines in their spectra. Since the redshifts 
of arcs Bl, B2, B3, and B4 are not known, an assumption 
needs to be made of whether all those arcs are images of one 
source or more sources. Based on very similar color, struc- 
ture and distance from the center of the cluster we suppose 
that arcs Bl, B2, B3, and B4 are images of one source and 
we call it system B. This conjecture is supported by the 
lensing model described below, since by assuming the exis- 
tence of two separate systems (B1-B2, B3-B4) our model 
predicts relatively bright multiple images which are not ob- 
served. We excluded that the feature D is a radial arc, de- 
spite its elongated morphology, since no tangential counter 



Fig. 1. The RCS0224-0002 cluster with labeled arcs. Color 
image composed from F814W and F606W WFPC2 HST 
images. The image is 40 arcsec across. 



images are visible and because its position and morphology 
makes this hypothesis unlikely. Our model suggests that 
feature C is a central demagnified image, which is clearly 
visible in Fig. |4] showing the F606W image after subtract- 
ing the two cD galaxies. There is also a very faint red arc, 
labeled E, which was not included in our analysis. 

Sinc e mass is known to follow light in galaxy clusters 
(see e.g- ISand et al-ll2002f ). the distribution of color selected 
cluster members is often used to model substructure of the 
underlying dark matter. Besides to the two brightest central 
galaxies (BCGs), there is no pubfic spectroscopic informa- 
tion available in the field, we then used the red sequence to 
identify likely cluster members. In Fig. [2] we show the color- 
magnitude diagram over the whole WFPC2 field, highlight- 
ing red sequence objects lying withing 15" from the clus- 
ter co re. Photometry was perf ormed using SExtractor soft- 
ware (jBertin fc Arnoutsl ll996V by detecting sources in the 
F814W band and measuring F606W - F814W colors with 
aperture of 1" diametefl The solid and dot-dashed lines 
represent ou r best fit to the red sequence and the best 
fit found by iBest et all (|2002D for the cluster MS1054 at 
z — 083 for the same filters, after applying a K-correction 
of 0.07 mag. Red sequence objects were defined as those 
within ±0.25 mag of the best fit line. 



4. X-ray emission 

The X-ray emission traces the hot gas trapped in the clus- 
ter potential well. The gas itself contributes about 15% 
to the total mass of the cluster and for relaxed systems 
traces closely the total mass density distribution. We over- 
lay the X-ray contours of RCS0224-0002 from the 100 
ksec Chandra observations in the 0.5-2 keV band onto 
the WFPC2 image in Fig. [5] The overall X-ray emission is 



The WFPC2 
to: ZPab 
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points were calculated accord- 
'2.5log{PHOTFLEM) - 21.1 - 
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Fig. 2. The color- magnitude diagram of RCS0224 with 
the WFPC2 F606W/F814W filters. The dots represent aU 
objects in the field. The squares represent the cluster red 
sequence (galaxies within 15 arcsec from the cluster cen- 
ter), the stars mark two central galaxies. The solid and 
dot-dashed lines are our best fit to the red sequence and 
the one of MS1054 at similar redshift. 




Fig. 3. The red sequence galaxies visible on the F814W 
filter WFPC2 HST image (objects that are marked here 
correspond to the squares in Fig [2l 



not symmetric, with a plume extending NW, and its peak 
shifted ~ 5 arc seconds north from the two central BCGs. 
To measure the X-ray temperature, we used and extraction 
region of 36.7 arcsec (or 265 kpc), which encompasses most 
of the X-ray emission by maximizing the signal-to-noise. 
The background subtracted, unfolded spectrum is shown 
in Fig. m We used Xspec v.12.3.0 Arnaud: (1996) to fit 
the data with a single te mperature Mekal model ( KaastTal 
Il992t iLiedahl et al.lll995^ and model the Galactic absorp- 
tion with tbabs (| Wilms et al.l 120001 ) . fixing the Galactic 
neutral Hydrogen columns de nsity to th e Galactic value ob- 
tained with radio data (Dicke v fc Lockm an 1990.). Since the 
signal-to-noise ratio in each energy bin is low, we used the 
C-statistics for the best fit model, over the energy range 0.6- 
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Fig. 4. RCS0224-0002 in the F606W filter with subtracted 
cD galaxies. The central radial feature C is clearly visible. 

8.0 keV (excluding low energy photons due to uncertainties 
of ACIS calibration) . We used 742 ± 35 total net counts in 
the fit (514 ± 23 in the soft 0.5-2 keV band) and found 
a best fit temperature of kT = 5.26'^\ q'^ keV (1-sigma 
error). The de-absorbed fiux within the extraction aper- 
ture, in the (0.5 - 2.0) keV band, is 1.84 x lO^^Vg cm^^ 
s~^ and the rest-frame X-ray luminosity Lx{0.5 — 2keV) — 
(0.38 ±0.02) X 10''^ erg s~^. The bolometric luminosity re- 
turned by the best fit model is Lbol = (1.28 ±0.06) x 10'*'^. 
With these values of X-ray luminosity and temperature, 
we note that RCS0224-0002, which is an optically selected 
cluster, lies on the Lx — T relation det ermined from large 
samples of X-ray selected clusters (e. g. iRosati et al.ll2003) 
We can use the measured cluster temperature to esti- 
mate the cluster mass assuming the hydrostatic equilibrium 
and isothermal distribution of the gas, with a polytropic in- 
dex 7 = 1. Using the standard /3- model for the gas density 
profile, Pgasir) — po/[^ + (r /re)^]^^^^, th e mass within the 
radius r can be written as (|Sarazinlll988f) : 



M(< r) ~ 1.11 X 



T{r) 



{r/7 



keV /i-iMpc 1 + (r/rc)' 



(1) 



A fit to the X-ray surface brightness profile with the cor- 
responding /3-model SB(r) oc [1 + (r/rc)^]^^^^"'^/^ yields a 
core radius = (253±72)kpc and (3 = 0.97±0.3. Therefore 
the mass within R200 = 0.4 Mpc is (1.7 ± 1.1) x lO^^M©. 



5. Model 

We constructed the mass model of RCS0224-0002 by fit- 
ting the position and shapes of the multiple image systems 
A, B and C. Based on the light distribution of most lu- 
minous red-sequence galaxies, our model consists of sev- 
eral mass components: two isothermal non-singular ellip- 
soids to reproduce global cluster properties (NIEl, NIE2); 
eight isothermal non-singular spheres fixed at the position 
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Fig. 5. The X-ray emission contours of RCS0224-0002 
(smoothed with a Gaussian with a — 5") over-plotted on 
the F606W WFPC2 HST image. 




Fig. 6. X-ray spectrum of RCS0224-0002 from lOOksec 
Chandra observations, with the best fit Mekal model, for 
kT = ?>.2&t\lj keV. 



of cluster members (NISI. .8) - refereed to as the substruc- 
ture; one non-singular ellipsoid, corresponding to the elon- 
gated object marked D in Fig. 1 (NIE3). In order to reduce 
the number of free parameters, we fixed the positions and 
the relative masses of the galaxy cluster clumps using the 
optical data available. In summary, we have 17 adjustable 
parameters in our model, including sources positions and 
unknown redshifts. All parameters are listed in Tables [1] 
and [21 The seven observed extended images are enough to 
constrain those 17 parameters due to the fact that we base 
our goodness of fit function not only on the position of the 
images but on the full information encoded in their shapes. 
Models including the radial feature D as a counter-image 
give the worst results, but as mentioned in Sect. [31 it is 
probably an foreground edge-on galaxy. Arc E was not used 
in the model since its redshift is unknown and it is too faint 
to provide any further constraint. We would like to empha- 



size that we do not assign any physical meaning to the two 
distinct smooth components (NIE1&NIE2), and we are in- 
terested in the properties of the overall, combined profile. 
We have also tried to fit the data with only one smooth 
component (NIEl) and the substructure, however in that 
case we were not able to fit the arcs system B accurately. 

5.1. Mass profiles 

Although the N-body simulations of dark matter halo 
formation suggest NFW profiles rather than isothermal 
ones, recent strong lensing studies do not exclude and 
i n some cases even prefer isothermal p rofile over NFW 
([Gavazzi et al.l 120031 : iHalkola et a^l2006^ . We model here 
all mass components as non-singular isothermal ellipsoid, 
a simple generalization of a non-singular isothermal sphere 
often used as a physical representation of a gravitationally 
relaxed system. The use of isothermal profiles has also the 
advantage of being computationally less demanding. The 
associated gravitational potential 4>, projected mass den- 
sity p, and deflection angle a are given by 



ai(a;i,a;2) 

a2ixi,X2) 

pixi,X2) 



xit: 1- X2t: — 

oxi 0x2 

- bqs\n[{^/j + s)^ + {1 - q^)xl]'^^ 

+ bqs\n[{l + q)s] , 

d(j) bq 
dxi J\ - qr2 



(2) 



dx2 

,2/„2 



bq 



a/1 - q'^Xi 



atanh 



q'is'+xj)+xl , 

1 b 

2 (s2 +xj+ x\lq^y. 



y/l - q'^X2 
tp + q^s 



(3) 

(4) 

(5) 
(6) 



where q is the ellipticity, s is the core radius, b is the scale 
factor and x = (xi,a;2) is the position in the image plane. 
Note that the fast method for calculating the deflection 
an gle of the soften ed non-singular ellipsoid might be found 
in iBarkanal (|1998[ ). 

5.2. Minimization mettiod 
5.2.1. Source plane minimization 

In order the get a first, approximated solution, we perform 
model fitting minimization on the source plane. This tech- 
nique is computationally very efficient, since there is no 
need to solve the inverse problem of the lensing equation 
and the deflection angle is only computed at the position 
of the images. We also assume that sources are small com- 
pared to the scale of variations of the lensing potential. If 
we have N images at positions {x{\ corresponding to one 
source, then we define the as 



Xsrc 



where 



' 6uJ pifS^ ^pidui + P 



Sui 



Mobs 
•^obs 



(7) 



(8) 
(9) 



J. Rzepecki et al.: Strong Lensing Analysis of the Cluster RCS0224-0002 at 2: = 0.77 



5 



In the previous equations, itobs,i is the source position (as 
predicted by the model) corresponding to the image a;obs,i, 
0(a;obs,i) is the lensing potential at i mage i, w is the cos- 
mological weight of the source (see e.g. lLombardi &: BertinI 
I1999D . and fii is the magnification matrix (inverse of the 
Jacobian matrix of the lens mapping) at the image i. 
Moreover, in Eq. ([7]) we introduced the covariance matrix 
of the position measurements S. For simplicity, in this pa- 
per we assume that the covariance matrix is diagonal and 
takes the form 



or 






^2 



(10) 



where ct^ is estimated to be ~ 0.05". In the definition of our 
[Eq. ([7])] we introduced also a "penalty" function P. This 
function, is used to bound some of the free parameters to 
certain intervals, and is chosen to have the functional form 



N 



-'down,?' 



ft)) 



P = P^atan(lO^(6d 

+ atan(10'^(pi - 6up,i)) + Vtt , 



(11) 



where, N is the number of bounded parameters in our 
model. Pi is the i-th bounded parameter, which is required 
to be in the range [6down,i, &up,j]- Note that the penalty 
function P behaves similarly to a "square potential well" , 
i.e. the sum of two Heaviside functions; however, the use 
of analytic functions ensures that P is differentiable and 
makes our minimization numerically stable. In order to 
effectively bound our parameters, we used a large number 
for the coefficient V. When there is more then one source, 
the same procedure is repeated for all sources and the 
resulting from Eq. ^ are added. The magnification 
matrix /i^ is included because HiSui ^ 6xi, so that Xsrc 
is an approximation of in the image plane. However, 
this also introduces a weight in the term, as images for 
which fiidui are small do not contribute significantly to the 
minimization process. It is possible to write an analytical 
expression for the source position that minimizes Xsi-c- 



Umod = A ^b;, (12) 

^=E^^^rV., (13) 

i 

b = ^ S:^^^ fj.iUohs,i ■ (14) 



5.3. Extended images 

The best fit model provided by Eq. ([7]) is used as start- 
ing point for the image plane analysis. This step is based 
on a new minimization, with a x^ composed of two 
terms: the so-called M odified Hausdorff Distance (MHD, 
iDubuisson fc Jainll994j) between the modeled and observed 
image sets and the "plain difference" between the same sets. 



The MHD between two sets A and B is defined as 
MHD = max(/iab,/iba) , 

" " aeA 



, , mm lla — 611 
LB I ^ aeA" " 



(15) 
(16) 

(17) 



In addition, the "plain difference" between the observed 
and modeled arcs is computed as follows. All pixels in each 
observed arc system, generically called O, are assigned a 
value of 1; other pixels are assigned a value of —1. The same 
procedure is applied to the corresponding modeled arcs (M) 
and the difference diff(0, Af) = \0 — M\ is calculated. In 
summary, the expression to minimize in the image plane is 



= MDH(i:>, M) + Lo diff (A/, D) + P 



(18) 



The factor lo was chosen to be ~ 0.1, since this value re- 
sulted in the fastest convergence. The penalty function P 
is used to bound some of the model parameters and it is 
defined in Sect. 15.2.11 By using two distance components, 
we ensure an efficient convergence of the minimization since 
when the modeled and the observed images start to overlap, 
the MHD becomes less sensitive to sma ll variations t hen the 
plain difference. The Powell algorithm (jPowellll 1964( 1 is used 
for all the minimization procedures. 



6. Results 

The best fit model (with MHD as defined by Eq. ([15]) equal 
to 30.3) is presented in Fig. [T] The values of corresponding 
parameters are given in Tables [T] and [2l The model repro- 
duces fairly well all the observed strong lensing features. 
The giant arc A include a counter-image 7" to the west of 
the BCGs (A3). The model also reproduce the quadrupole 
system B (B1,..B4). The central feature C is also predicted 
fairly close to the observed one, although with different 
morphology. None of the models we analyzed could repro- 
duce the radial feature D, which suggests that it is probably 
a foreground edge-on galaxy. In addition, inclusion of D to 
the lens model (NIE3) significantly improved our fits and 
allowed us to "break" the arcs system B into two arcs Bl 
and B3. The best fit redshift of the source for the system B 
is 2.65 ± 0.08; a spectroscopic redshift of these blue arcs, as 
well as object D, would provide a strong validation of our 
lensing model and could also be used to better constrain the 
mass distribution. Estimates of the statistical errors are dis- 
cussed in the following section. Figure [5] and Tab. [3] show 
the results of some tests performed to assess how well the 
best fit model is able to reproduce the morphology of the 
multiple image systems A and B. For this purpose, we ray- 
trace a given image for each system (A2 and Bl, marked 
with green boxes in Fig. [H]) into the source plane by using 
its HST color image. This gives us the reconstructed source 
image. We then ray-trace back all the pixels from the source 
plane into the image plane, thus finding all counter-images 
of the given image. These reconstructed counter-images are 
finally compared with the observed ones (Al,3 and B2,3,4). 
In general, we find a good agreement, especially the knots in 
the Al arc are very well reconstructed. The overall shapes of 
all the arcs in the system B are also accurately predicted. 
The mass of the cluster within -R200 = 0.4Mpc obtained 
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from the model is 1.9 ± 0.1 x 10^^ M© and its distribution 
is shown in Fig. [S) This is in a good agreement with an mass 
derived above from the X-ray temperature. Since we do not 
know all the cluster member galaxies, we cannot reliably es- 
timate the mass-to-light ratio of the whole cluster. For the 
substructure (the mass associated with the luminous cluster 
component - NISI. .8), we find an average mass-to-light ra- 
tio M/LB,vega ~ 3.6 Mq/Lo^b- Wc Converted the observed 
F814W filter flux to the rest frame B filter flux, by calcu- 
lating a k-correctiq n for a template elliptical galaxy from 
iKinnev et al.l (119961 ). The center of the mass of the best fit 
model follows the light distribution. NIEl is found to be a 
diffuse (core radius « 15 ") mass component close to the 
peak of the X-ray emission. The latter is shifted « 5" from 
the NIE2 component, which corresponds to the center of 
the potential well and the position of the BCGs. This may 
indicate the presence of a merger. The radial average pro- 
file of the best fit surface mass density is shown in Fig. [101 
This can be well approximated by a power law profile with a 



slope 7 = 0.74 



+0.03 
-0.04' 



which is closer to the isothermal profile 



(7 = 1) than results obtained in other clusters. For exam- 
ple, the analysis of the cluster J1004-4112 yielded 7 0.5 
(ISharon et a l.l l2005D and 0.3 < -y <0^ ()WiUiams fc Sahal 
120041 ) ■ whereas iBroadhurst et all (|2005aD found 7 = 0.5 in 
A1689 using a large number of identified multiple images. 
Note that the flat core of the mass profile we have found, 
being a result of a hight value of the Vc of the NIEl compo- 
nent, is well constrained by the position of the central arc 
C. The change of the Vc by 50% causes the shift in the C 
arc position of « 1 arc sec. 

By approximating the mass density distribution with 
NFW-like profile of the form 



7. Error analysis 

Our method involves the minimization of the MHD whose 
expression (Eq.18) is not a formal and includes a number 
of penalty functions (weights) to limit the range of some pa- 
rameters. As a result, it is difficult to obtain reliable errors 
on the best fit parameters. In the presence of many pa- 
rameters , the Mont e Carlo Markov Chain (MCMC, see for 
example 'Neall 1 1 9931 ) method is an efficient way to estimate 
the likelihood associated to our best fit model. MCMC is 
used as a third step of our minimization process by recon- 
structing the probability distribution function of our model 
parameters. We start the const ruction of Markov cha in us- 
ing the Metropolis algorithm (Metropo lis et al] |1953l from 
the best fit solution of the MHD minimization. We use a 
number of chains randomly distributed around the best fit 
point. The resulting chain being the composition of all those 
partial chains provides an approximate probability distri- 
bution function for our parameters, from which we estimate 
the confidence levels shown in the Fig.[T3 Also by randomly 
probing the parameters space, the MCMC algorithm helps 
to fine tune our best fit parameters returned by the previ- 
ous step of minimization. Most of the parameters are well 
constrained (within 10 - 20 percent). The unknown red- 
shift of the arc system B appears to be well constrained, 
zb = 2.65±0.08. The mass to light ratio of the substructure 
is however poorly constrained to be 3.6l^'g. 

We estimated the errors of a single power law and NFW- 
like profile parameters by drawing a random sample of mod- 
els from our Markov Chain, and then fitting a single power 
law and NFW-like profile to that sample. The resulting 
error estimates are presented in Fig. [T31 This shows that 
isothermal and NFW profiles are excluded with 99% confi- 
dence level. 



(r/rc)/5(l + r/rc)(i-« 



we find a slope 0.69„q'"3, flatter then the canonical 
NFW model (/3 = 1), however in good agreemen t with other 
studie s which obtained P < 1. For example, ISand et all 
(|2002D finds P = 0.35 for the galaxy cluster MS1237-23, 
and (3 < 0.57 (at 99% confide nce level) from the analysis of 
a large sample of clusters (,Sand et al.ir2004' ). 

In addition, we have tried to fit a model based the uni- 
versal NFW profile rather than NIE. The result, presented 
in the Fig. [11] shows that an NFW model performs signif- 
icantly worse then the NIE one. The arcs Al and A2 are 
reproduced fairly well, but the counter image A3 is found 
much too far from the cluster center. In addition, in the 
NFW model feature B4 is split into two arcs (the second 
of which is not observed) and the reproduced arc B2 is 
shifted with respect to the observed one. This is reflected 
by the value of MHD, which is ten times bigger then the 
corresponding value for the best-fit NIE model. We note, 
however, that this bad performance might be due to the ap- 
proximated NFW elliptical model used in our code, where 
the ellipticity is achieved by perturbing the potential of the 
spherical NFW profile instead of its density. This approxi- 
mation holds for potentials close to spherical, and therefore 
we need to impose additional restrictions on the ellipticity 
of the NFW components. 



8. Conclusions 

We have performed a strong lensing analysis of the cluster 
RCS0224-0002 using HST/WFPC2 images in F814W and 
F606W bands. We used two arc systems: a red giant tan- 
gential arc 14" from the center, with measured redshift of 
4.87, for which we identified an inner counter image, and a 
system of blue arcs at smaller radii with no spectroscopic 
information. 

We have modeled the mass distribution with with three 
mass components: isothermal spheres associated with the 
most luminous cluster members to model the substruc- 
ture, and two isothermal ellipsoids to model the underlying 
smooth mass component. Since spectroscopic information 
is available in the literature only for two cD galaxies, we 
identified likely member galaxies in the cluster core from 
the red sequence, which is clearly detected in the F606W- 
F814W color distribution. To infer the mass distribution 
from the position and shapes of the strong lensing features 
we used a three-step approach: i) minimization of the size of 
the two sources on the source plane, ii) minimization of the 
difference between the observed and modeled arcs on the 
image plane, based on the Modified Hausdorff Distance, 
and iii) a refined estimate of the best fit parameters and 
errors analysis with the Monte Carlo Markov Chain. The 
resulting mass density reproduces all the strong features 
fairly well. The redshift of the blue arc system is predicted 
to be 2.65 ±0.08. 
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We find that the substructure made of nine isother- 
mal components centered on the brightest cluster members, 
with M/LB.t,ega ~ 3.6Mq/Lq^b is crucial to exactly repro- 
duce the shapes and positions of all the arcs. 

By fitting a single power-law or NFW-like halo to the ra- 
dial average mass density distribution we have found that 
both profiles are far from canonical isothermal and stan- 
dard NFW: we have found the power-law parameter 7 to 
be 0.74^Q Q4 (7 = 1 for an isothermal profile) and steepness 

parameter for NFW-hke profile /3 to be 0.69^^3 (/^ = 1 
for a NFW profile), with the upper boundary very well 
constrained. Both those values are consistent with the re- 
sults obtained by stu dying the strong l ensing properties 
of other clusters (see ISand et all l2002l 120041 ). The best 
fit NIS has dv = 925km/s and rc = llkpc; the best fit 
NFW has i?200 = 0.4Mpc and concentration parameter 
c = 3.4tQ;5, similarly to other massive c lusters (c w 4 
for a z = 0.18 cluster iHal kola e t alll2006l . c w 5 for z = 
0.68 cluster IWilliams fc Sah a 200^. However, a wide range 
of concentr ation parameter s are found (e.g. for c > 10 
see iBroadhurst et al.ll2005bl ). We have measured the total 
mass of the cluster within i?200 to be 1.9 ± 0.1 x IQ^'^Mq 
and its main component may be well described by a two 
NISs with a a^i = 94512° km/s, a r^i = 112+}i^kpc, a 
(^v2 = 702^28 km/s, and a — 12^2 kpc. The mass 
of RCS0224-0002 derived from the lensing model is in a 
very good agreement with the one obtained from the X- 
ray temperature measured with deep Chandra observations 
(M2oo = (1.7±1.1) X IQI^Mq). 

This analysis shows that even with a limited number 
of identified multiple images we could constrain the mass 
distribution fairly accurately. This was possible, in the case 
of RCS0224-0002, because the two arcs systems are at very 
different angular diameter distances and probe significant 
fraction (« 20% for the arcs system A, and ~ 60% for the 
system B) of the Einstein rings. Further spectroscopic ob- 
servations of the system B, as well as cluster members, will 
allow a very robust constraint of the mass density profile 
of the inner core of this cluster and its substructure. 

Acknowledgements. We would like to thank Matthias Bartelmann for 
very useful discussions and comments, specifically on the usage of the 
MCMC method to estimate the errors of our model. 
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Table 1. Parameters defining our model (see equation [^Hl]) after minimization. Parameters in parenthesis were allowed 
to change during minimization 





NIEl 


NIE2 


NISI 


NIS2 


NIS3 


NIS4 


NIE3 


NIS5 


NIS6 


NIS7 


NIS8 


Xi 


(16.834) 


19.413 


18.039 


20.578 


24.799 


23.494 


17.621 


22.389 


23.614 


25.296 


14.097 


X2 


(18.502) 


20.307 


20.834 


20.147 


7.336 


15.330 


10.253 


12.612 


12.005 


32.969 


16.076 


z 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


0.782 


b 


(19.196) 


(10.086) 


(0.116)* 


(0.191)* 


(0.027)* 


(0.032)* 


(0.417) 


(0.081)* 


(0.075)* 


(0.087)* 


(0.037)* 


q 


(0.396) 


(0.597) 










0.3 













(2.994) 


(1.409) 










0.873 










s 


(15.539) 


(1.778) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 



Xi,X2 : central position in arc seconds in the coordinate system of the Fig. |21 z : redshift, b : scale factor in arc seconds, q : 

ellipticity, 9 : position angle in radians, s : core radius in arc seconds 
* - for the substructure the M/L ratio has been used as the variable for the minimization 



Table 2. Parameters defining sources after minimization. Parameters in parenthesis were allowed to change during 
minimization 





SOURCEl 


S0URCE2 


Ul 


(15.979) 


(18.892) 




(19.204) 


(19.412) 


z 


4.878 


(2.648) 



Ul, U2 source position in arc seconds, z : redshift 



Table 3. Images reproduced under image plane - source plane - image plane mapping 
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[ore sec] 

Fig. 7. Images reproduced by our best fit model over-plotted on the combined F606W/F814W WFPC2 HST image. The 
closed lines show the critical curves and caustics for a source at z — 4.87. The center of the image is at RA 02:24:34.218, 
Dec -00:02:31.64. 




Fig. 8. Result of image plane ~ source plane - image plane mapping. Panel to the left shows the arcs (marked by boxes) 
used to reproduce the arc systems. Middle and right panels show the arc systems as reproduced by the best fit mass 
model (the box marks the original image). 
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Fig. 9. Mass density produced by our best fit model. The closed lines are the critical curves for a source at z = 4.87. The 
crosses (+) mark the positions of our model components. The big cross (X) gives the position of the peak of the X-ray 
emission. The center of the image is at RA: 02:24:34.218 Dec: -00:02:31.64, the orientation as in Fig. [5] 
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Fig. 10. Radial average profile of the surface mass density of our best fit model (solid line) versus power law profile with 
7 = 0.74 (dotted line) and NFW-hke profile with /3 = 0.69 (dashed line). 
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Fig. 11. Images reproduced by our best fit NFW model over-plotted on the combined F606W/F814W WFPC2 HST 
image. The closed lines are the critical curves for a source at z = 4.87. The center of the image is RA: 02:24:34.218 Dec: 
-00:02:31.64 
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Fig. 12. MCMC error estimates. The contours correspond to 68%, 90% and 99% confidence levels. Marks on vertical 
and horizontal axis give the same confidence levels for ID projected variables. The cross marks the position of the best 
fit point. 
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Fig. 13. MCMC error estimates of the parameters of the single power law and NFW-like profiles fit. The contours 
correspond to 68%, 90% and 99% confidence levels. Marks on vertical and horizontal axis give the same confidence levels 
for ID projected variables. The cross marks the position of the best fit point. 



